Strong Decays of Light Vector Mesons 
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The vector meson strong decays p — > niv, <j> — > KK, and K* — > nK are studied within a covariant 
approach based on the ladder-rainbow truncation of the QCD Dyson-Schwinger equation for the 
quark propagator and the Bethe-Salpeter equation for the mesons. The model preserves the one-loop 
behavior of QCD in the ultraviolet, has two infrared parameters, and implements quark confinement 
and dynamical chiral symmetry breaking. The 3-point decay amplitudes are described in impulse 
approximation. The Bethe-Salpeter study motivates a method for estimating the masses for heavier 
mesons within this model without continuing the propagators into the complex plane. We test the 
accuracy via the p, <f> and K* masses and then produce estimates of the model results for m ai and 
mt 1 as well as the proposed exotic vector 7Ti(1400). 
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I. INTRODUCTION 

The study of light-quark pseudoscalar and vector 
mesons is an important tool for understanding how QCD 
works in the non-perturbative regime. The pseudoscalars 
are important because they are the lightest observed 
hadrons and are the Goldstone bosons associated with 
dynamical chiral symmetry breaking. The ground state 
vector mesons are important because, as the lowest spin 
excitations of the pseudoscalars, they relate closely to 
hadronic qq modes that are electromagnetically excited. 
The decay of vector mesons to a pair of pseudoscalars 
proceeds via a P-wave interaction; this probes a differ- 
ent aspect of the meson Bethe-Salpeter [BS] amplitudes 
than does the electroweak decay constant which is essen- 
tially the projection of the relative wavefunction onto the 
origin of separation. 

To calculate these coupling constants, we use an ap- 
proach based on the Dyson-Schwinger equations [DSEs] , 
which form an excellent tool to study nonperturbative as- 
pects of hadron properties in QCD jj]]. The approach is 
consistent with quark and gluon confinement [jR 0, ^ |J, 
generates dynamical chiral symmetry breaking |jfp| , and 
is Poincare invariant. It is straightforward to imple- 
ment the correct one-loop renormalization group behav- 
ior of QCD fi, and obtain agreement with perturba- 
tion theory in the perturbative region. Provided that 
the relevant Ward-Takahashi identities [WTIs] are pre- 
served in the truncation of the DSEs, the corresponding 
currents are conserved. Axial current conservation in- 
duces the Goldstone nature of the pions and kaons ||; 
electromagnetic current conservation produces the cor- 
rect electric charge of the mesons without fine-tuning. 
These properties are implemented here within the rain- 
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bow truncation of the DSE for the dressed quark prop- 
agators together with the ladder approximation for the 
Bethe-Salpeter equation [BSE] for meson bound states. 
The model we use has two infrared parameters. Previ- 
ous work has shown this model to provide an efficient 
description of the light-quark pseudoscalar and vector 
mesons ]?], |[|. Furthermore, in impulse approximation, 
the elastic charge form factors of the pseudoscalars 
and the electroweak transition form factors of the pseu- 
doscalars and vectors p|, [l2j are in excellent agreement 
with data. This suggests that the strong decays of the 
vector mesons should be well-described in impulse ap- 
proximation without parameter adjustment. 

The Euclidean metric that we employ facilitates the 
modeling of the gluon 2-point function or more generally 
the effective quark-quark interaction; it also is the metric 
within which practical DSE solutions are available in the 
literature. A complicating element is that solution of the 
BSE for meson bound states requires an analytic contin- 
uation in the meson momentum to reach the mass-shell. 
The resulting quark p 2 in the Bethe-Salpeter integral 
covers a certain domain in the complex plane that differs 
from the real spacelike axis by an amount that increases 
with meson mass. For low mass mesons such as tt and K, 
the required continuation of DSE solutions for the quark 
propagator is not difficult. For the ground state vector 
mesons, the difficulties are much greater but have been 
overcome directly ||. For meson masses greater than 
about 1 GeV 2 the numerical requirements of this proce- 
dure become burdensome. Furthermore, complex plane 
singularities of rainbow DSE solutions for propagators 
limit this direct approach to mesons with masses below 
about 1.2 GeV JT|. 

Here we test the reliability of an approximation that 
uses a Taylor expansion to make the continuation away 
from the real spacelike p 2 axis. Since information about 
quark propagators is then needed only on the real space- 
like p 2 axis, this approximation is considerably easier to 
implement. It is similar in spirit to an approximation 
employed some time ago [|3 . If the behavior for real mo- 
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menta dominates the physics, low orders of a Taylor ex- 
pansion should be effective. This work is an exploration 
of that conjecture since we can test against the direct 
complex plane results || for BSE solutions for masses 
of p, K* and <p within the same model. Independently 
of how the behavior of the propagators in the complex 
plane is obtained, the calculation of meson form factors 
or decays in impulse approximation requires that the rel- 
ative quark momentum of each vertex amplitude be also 
continued to the complex plane. In the interests of con- 
sistency, we employ a Taylor expansion method for both 
vertex amplitudes and propagators in the calculation of 
decays. 

The results from the Taylor expansion method are used 
to motivate a method for estimating the masses and BSE 
solutions for heavier mesons within this model without 
continuing the propagators into the complex plane. We 
test the accuracy of this real-axis approximation for the 
ground state pseudoscalar and vector mesons, and then 
produce estimates of the model results for m ai and 
as well as the proposed exotic 1 h vector 7r 1 (1400). 

In Sec. II we review the formulation that underlies a 
description of vector meson strong decays within a mod- 
eling of QCD through the DSEs and BSEs. We discuss 
the diagrams needed for a calculation in the impulse ap- 
proximation for p — > 7T7T, (f> — > KK, and K* — > Kit. In 
Sec. Ill we discuss the details of the model and we in- 
clude a comparison with recent lattice-QCD results for 
the quark propagator. In Sec. IV we present our nu- 
merical results for masses and decays, and discuss the 
performance of the calculational technique. In Sec. V 
the technique for estimating the masses of heavier states 
is evaluated. A concluding discussion is given in Sec. VI. 



II. 



VECTOR MESON STRONG DECAYS 



We use the Euclidean metric where {7^,7^} = 28^ 



7, 



t - 



7^ and a ■ b — ~^2 i=1 oabi. The invariant amplitude 
for the coupling of a vector state with helicity A and mo- 
mentum P1+P2 to a pair of pseudoscalars (or scalars) hav- 
ing momenta pi/2 has the form M. A = A^e^ where e A is 
the polarization vector of the vector state, and the vertex 
which describes the coupling of the Lorentz component 
of the vector state to the pair of pseudoscalars has the 
form A M = g (p\ — P2)p, where g is the coupling constant. 
In the case that the vector state and the pseudoscalar 
states are qq bound states with flavor labels ab, ca and 
be respectively, the impulse approximation for A^ b,c can 
be expressed as 



vector, and P = (pi — p 2 )/2. We denote the properly 
normalized BS amplitude at a vertex having an outgoing 
quark of flavor /1 and momentum p + and an incoming 
quark of flavor f% and momentum p_ by Tf 1 ^ 2 {p +1 p^). 
The notation J, = J d A k/(2n) A stands for a transla- 
tionally invariant regularization of the integral, with A 
being the regularization mass-scale. The regularization 
can be removed at the end of all calculations, by taking 
the limit A — * 00. It is understood that the same regular- 
ization is applied to the calculated dressed quark propa- 
gators 5(g) and BS amplitudes appearing in Eq. and 
that the quark propagators are renormalized at a conve- 
nient spacelike momentum scale before taking the limit 
A — > 00. 

For mesons with a unique quark flavor content, the 
physical vertices are given by the A£ fc < C (P; Q). For states 

like \p° >= (\uu > — \dd>)/y/2, linear combinations are 
required. We have, for example, 



A" 



A 



K* 4 



*K+K- 



A 



= ±[Ap d (P;Q) 

= A^"(F;Q), 
= ±A^(P;Q), 



= K 



d (P;Q) 



A dd >"(-P;Q)],(2) 

(3) 
(4) 
(5) 



The vertices for modes that differ from these by electric 
charge are obtained by quark flavor exchange. For exam- 
ple under isospin exchange of u and d flavors, the last two 
modes become K*° — > K°tt° and K*° — ► K + ir~ respec- 
tively. In this work we use isospin symmetry in which u- 
and G?-quarks have identical strong interactions and prop- 
erties except electric charge. Under these circumstances, 
A"'" = A ds,d = A™ s ' d and there is only one independent 
K* — > Ktt vertex. 

The physical coupling constants for each v — > pp decay 
are identified from the vertices at the mass shell according 
to 



A^(P;Q2 = _ m 2) = 2 PT 9v ^ 



(6) 



where P T is the component of P perpendicular to the 
vector meson momentum Q. The decay width of a vector 
of mass m v is given by 



A=l 



(7) 



where p is the invariant phase space factor, the invari- 
ant amplitude is A4 X = A p and e A is the polarization 
vector. The explicit form is 



A^ C (P;Q) = tr sc / 5 c (g) T%(q, q+) S a (q+) 

Jk 

xr^(g+,g_)5 b (g_)f^(g_,g), (1) 

where q = k + P/2 and q± — k — P/2 ± Q/2, and tr sc de- 
notes a trace over color and Dirac spin. For convenience 
we have introduced Q = p\ + p 2 as the momentum of the 
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A ' {m v ,m x ,m 2 ) , 



(8) 



where 



X(a,b,c) = a z + b z + c z - 2ab~ 2ac~ 2bc , (9) 

and m p i, and m P 2 are the masses of the pseudoscalar 
decay products. 
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A. Dyson— Schwinger Equations 

The dressed quark propagator and the meson BS am- 
plitudes are solutions of their respective DSEs, namely 

+ Z x j g 2 D tlv (jp-q)^S{q)Yi{q,p) , (10) 

and 

T ah (p + ,p_) = ^ K{p,q-Q)S a { q+ )T a \q + ,q^)S\q^), 

Jq 

(11) 

where D^^k) is the renormalized dressed-gluon propaga- 
tor, T]j{q,p) is the renormalized dressed quark-gluon ver- 
tex, and K is the renormalized qq scattering kernel that 
is irreducible with respect to a pair of qq lines. The quark 
momenta are q± ; the meson momentum is P — q+ — q- 
and satisfies P 2 = — m 2 . 

The solution of Eq. (jl(^) is renormalized according to 
S(p)^ 1 = ij-p + m(p) at a sufficiently large spacelike [X 2 , 
with mUi) the renormalized quark mass at the scale \i. 
In Eq. (|ll|), S, r* , and m(/z) depend on the quark flavor, 
although we have not indicated this explicitly. The renor- 
malization constants Z 2 and Z± depend on the renormal- 
ization point and the regularization mass-scale, but not 
on flavor: in our analysis we employ a flavor-independent 
renormalization scheme. 

To implement Eq. (|ll|) while preserving the constraint 
P = q+ — q- , it is necessary to specify how the total mo- 
mentum is partitioned between the quark and the anti- 
quark. The general choice is q+ = q + rjP and g_ = g— 
(1 — rj)P where rj is the momentum partitioning param- 
eter. The choice of rj is equivalent to a choice of relative 
momentum q; physical observables should not depend on 
the choice. This provides us with a convenient check on 
numerical methods. 

The meson BS amplitude r a& is normalized according 
to the canonical normalization condition 

2P M = tr sc -^-{ J r ai (q',q) S a (q + )T a \q,q') S b (q^) 
+ [ A x ah Ck',~k)K(k,q;P) X a \q,q')} , (12) 

Jk,q ' 

at the mass shell P 2 = Q 2 = — m 2 , with q = q+ i]Q, 
q' = q — (1 — f])Q, and similarly for k and k'. Thus the 
derivative in Eq. ( |l2| ) acts only upon the propagators S 
and the kernel K. For vector mesons, it is understood 
that one must contract and average over the Lorentz in- 
dices of the (transverse) BS amplitudes due to the three 
independent polarizations. 

For pseudoscalar bound states the BS amplitude is 
commonly decomposed into Q 

r(<z+,g_) =7 5 [iE+ f F+ iG + a^q^H] , (13) 



where the 4 independent covariants have been con- 
structed from the momentum basis q and P rather than 
from g + and g_. Hence the invariant amplitudes E, F, 
G and H are Lorentz scalar functions f(q 2 , q ■ P;rj) that 
have an explicit dependence upon the momentum parti- 
tioning parameter r\. The dependence of the amplitudes 
upon q ■ P can be conveniently represented by the follow- 
ing expansion based on Chebyshev polynomials 

oo 

f(q 2 ,q-P) = X>i(cOS0)(gP)7i(« 2 ), (14) 

i=0 

where cos 8 = q ■ P/(qP). For charge-parity eigenstates 
such as the pion, each amplitude E, F, G, and H will 
have a well-defined parity in the variable q ■ P if one 
chooses rj = 1/2. In this case, these amplitudes are either 
entirely even (E, F, and H ) or odd (G) in q ■ P, and only 
the even (E, F, and H) or odd (G) Chebyshev moments 
fi are needed for a complete description. 

Since a massive vector meson bound state is trans- 
verse, the BS amplitude requires eight covariants for its 
representation. We choose the transverse projection of 
the form 

r^fe+.g-) = i f ,v 1 + q ^v 2 + q^Pv 3 

+ (Tfj, u P v Vj + q^CFapqaPp V S . (15) 

This form is a variation of that used in Ref. |J that is 
simpler and easier to use in many respects. The invariant 
amplitudes Vi are Lorentz scalar functions of q 2 and q ■ P 
and, for charge eigenstates, they are either odd or even 
in q ■ P. For the 1 rho meson, V3 and V§ are odd, 
the other amplitudes are even. To reverse the charge 
parity, one simply reverses this odd-even property of the 
amplitudes. 

III. LADDER-RAINBOW MODEL 

We employ the model that has been developed recently 
for an efficient description of the masses and decay con- 
stants of the light pseudoscalar and vector mesons [Q, |) . 
This consists of the rainbow truncation of the DSE for the 
quark propagator and the ladder truncation of the BSE 
for the pion and kaon amplitudes. The required effective 
qq interaction is constrained by perturbative QCD in the 
ultraviolet and has a phenomenological infrared behav- 
ior. In particular, the rainbow truncation of the quark 
DSE, Eq. ©, is 

Z ig 2 D^(k)Tl(q,p) - g(fc2 )jD f-(fc) 7 ^ (16) 

where D^ c (k = p — q) is the free gluon propagator in 
Landau gauge. The consistent ladder truncation of the 
BSE, Eq. (|n|), is 

K(p+, q + ;P) - -G(k 2 ) D^ c (k)^ ® , (17) 
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TABLE I: The pseudoscalar observables that define the 
present ladder-rainbow DSE-BSE model, adapted from 
Refs. g|. 





experiment |l5|] 


calculated 




(estimates) 


(t fitted) 


m ii=lGeV 


5 - 10 MeV 


5.5 MeV 


m /j=lGcV 


100 - 300 MeV 


125 MeV 


- (qq)l 


(0.236 GeV) a 


(0.24V) 3 


m w 


0.1385 GeV 


0.138 f 


u 


0.131 GeV 


0.131 1 " 


m K 


0.496 GeV 


0.497+ 


fx 


0.160 GeV 


0.155 



where k = p — q. These two truncations are consistent in 
the sense that the combination produces vector and axial- 
vector vertices satisfying the respective WTIs. In the ax- 
ial case, this ensures that in the chiral limit the ground 
state pseudoscalar mesons are the massless Goldstone 
bosons associated with chiral symmetry breaking J?], |). 
In the vector case, this ensures electromagnetic current 
conservation. 

The model is completely specified once a form is cho- 
sen for the "effective coupling" G(k 2 ). The ultraviolet be- 
havior is chosen to be that of the QCD running coupling 
a(k 2 ); the ladder-rainbow truncation then generates the 
correct perturbative QCD structure of the DSE-BSE sys- 
tem of equations. The phcnomenological infrared form of 
G(k 2 ) is chosen so that the DSE kernel contains sufficient 
infrared enhancement to produce an empirically accept- 
able amount of dynamical chiral symmetry breaking as 
represented by the chiral condensate . 

We employ the Ansatz found to be successful in earlier 
work J?], |J 



G(k 2 ) 



4tt 2 D k 2 



4n 2 lm T(k 2 



I In 



r+(l + P/A| CD )' 



(18) 



with 7„ 



12 



and T(s) = (1 - exp( 3 -^))/ S . The 



first term implements the strong infrared enhancement 
in the region < k 2 < 1 GeV 2 required for sufficient dy- 
namical chiral symmetry breaking. The second term 
serves to preserve the one-loop renormalization group 
behavior of QCD. We use m t = 0.5 GeV, r = e 2 — 1, 
Nf = 4, and we take A QC d = 0.234 GeV. The renormal- 
ization scale is chosen to be ijl = 19 GeV which is well into 
the domain where one-loop perturbative behavior is ap- 
propriate (7[ The remaining parameters, u> = 0.4 GeV 
and D = 0.93 GeV 2 along with the quark masses, are 
fitted to give a good description of (qq), m n /K and /„■. 
The subsequent values for fx and the masses and decay 
constants of the vector mesons p, (j>, K* are found to be 
within 10% of the experimental data |}|, see Tables | and 
III. A detailed analysis of the relationship between QCD 



and this Landau gauge, rainbow-ladder truncation of the 
DSEs with renormalization group improvement, can be 
found in the originating work JtJ . 
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In Figs. and | we compare the DSE model H propa- 
gator amplitudes defined by S(p) — Z(p 2 )[if> + M(p 2 )]~ 1 
with the most recent results in lattice QCD using stag- 
gered fermions in Landau gauge ]T7| , [i~8] |. These simu- 
lations were done with the Asqtad improved staggered 
quark action, which has lattice errors of order C(a 4 ) and 
0(a 2 g 2 ). In previous such comparisons |ll| ||(| the lat- 
tice data was less reliable. Fig. [j] shows both M{p) and 
Z(p) obtained with a bare lattice mass of ma = 0.036 in 
lattice units, which corresponds to a bare mass of 57 MeV 
in physical units. The DSE calculations use a current 
mass value of 75 MeV at fi = 1 GeV to match the lattice 
mass function around 3 GeV; this current mass is about 
0.6 m s . There is agreement in the qualitative infrared 
structure of the mass function particularly in the way the 
infrared enhancement sets in. Since the lattice simulation 
produces the regulated but un-renormalized propagator, 
the scale of the field renormalization function Z is ar- 
bitrary and only the shape is a meaningful comparison. 
For this reason, we have rescaled the lattice data for Z 
so that they match the DSE solution above 3 GeV. The 
ladder-rainbow DSE model typically produces a Z that 
saturates much slower than does the lattice Z] this may 
signal a deficiency of the bare gluon-quark vertex. 

The lattice work |l7], [l8| also produced a linear extrap- 
olation to the chiral limit mass function Mq(p) and we 
compare this to the DSE result in Fig ||. In principle 
the leading UV behavior of this mass function is propor- 
tional to the chiral condensate. The extraction of this 
from earlier lattice data has been attempted but there 
are uncertainties still to be resolved pl| . We note here 
that above 1 GeV there is excellent agreement between 
the present DSE model and the lattice results. 

Recent reviews put this model in a wider per- 

spective. These reviews include a compilation of results 
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FIG. 2: The chiral limit DSE mass function compared to the 
lattice chiral extrapolation |17J, n8| 



for both meson and baryon physics with similar models, 
an analysis how quark confinement is manifest in solu- 
tions of the DSEs, and both finite temperature and finite 
density extensions. The question of the accuracy of the 
ladder-rainbow truncation has also received some atten- 
tion; it was found to be particularly suitable for the flavor 
octet pseudoscalar mesons since the next-order contribu- 
tions in a quark-gluon skeleton graph expansion, have a 
significant amount of cancellation between repulsive and 
attractive corrections [E3| . 



IV. RESULTS 

A. Method of Calculation 

In order to identify the meson mass m from the BSE, 
Eq. (fTl|), in Euclidean metric, it is convenient to first 
introduce the linear eigenvalue A(P 2 ) of the kernel, con- 
tinue the equation in P 2 to the timelike region, and then 
find m such that A(— m 2 ) = 1. The BSE kernel involves 
terms linear in the 4-vector P and they are given the 
continuation P — > [im, 0). Thus with the quark propaga- 
tors written as S(p) = —ip'(Jv(p 2 )+ cs(p 2 ), the kernel of 
the BSE involves the amplitudes a a (q 2 _), where a = V, S 

and q\ = q 2 —m 2 /4 zkiqm z. Here q = yfq 2 > with 
q 2 being the (spacelike) integration variable, — 1 < z < 1 
is a direction cosine, and we have adopted equal mo- 
mentum partitioning rj = 1/2 for convenience. The am- 
plitudes cr a (q±) are in principle required to be known 
in a parabolic domain of the complex q± plane that 
includes the positive real axis and that extends sym- 
metrically in the imaginary direction and in the nega- 
tive (timelike) direction by amounts that grow with m. 
Previous work within the present approach and model 
has proceeded by use of the quark DSE to make the re- 
quired analytic continuation of these propagator ampli- 
tudes. For m > 1 GeV, the difficulty of the necessary 



TABLE II: Masses and decay constants (in GeV) for n and 
K obtained from the Taylor expansion treatment of q\ in 
quark propagator amplitudes. The effect of truncation of the 
4 Dirac covariants to just the canonical 75 covariant is indi- 
cated. Comparison is made with results from use of quark 
amplitudes along 5R(</±) only, and with the model-exact re- 
sults. 





m n 


u 


m K 


Ik 


1 st order Taylor 










All 4 ampls 




0.117 


0.111 


0.420 


0.134 


E ampl only 




0.106 


0.086 


0.383 


0.103 


2 nd order Taylor 










All 4 ampls 




0.138 


0.131 


0.498 


0.157 


E ampl only 




0.121 


0.098 


0.436 


0.116 


3 rd order Taylor 










All 4 ampls 




0.138 


0.131 


0.496 


0.153 


E ampl only 




0.121 


0.098 


0.435 


0.114 


Real axis only 




0.123 


0.099 


0.440 


0.127 


Model exact [9 




0.138 


0.131 


0.497 


0.155 


Experiment 


0.1385 


0.131 


0.496 


0.160 



numerical methods can outweigh the benefits of solution, 
and certainly accuracy becomes problematic. For this 
reason, we explore an approximate method. 

The eigenvalue A(P 2 ) should be real to obtain a real 
mass (no meson decay mechanisms are present in the lad- 
der BSE kernel) and thus one can speculate that the sub- 
domain consisting of the positive real axis for q\ might 
dominate the physics. To explore this systematically, we 
use a Taylor expansion of a a (q±) to reach the argument 
values off the positive real axis. Thus the only infor- 
mation needed about the quark propagator amplitudes 
from solution of the DSE are their values and derivatives 
on the Euclidean momentum domain. Within the BSE 
kernel we thus employ 

a a {q±) = a a (q%)+A 2 ± a' a (q 2 E ) + --- , (19) 

where q 2 E = max[0, ^(q 2 .)] and q E + = q±. For prac- 
tical reasons we are interested only in low orders of this 
expansion. 

B. Masses and Decay Constants 

In Table || we show the results for the masses and 
decay constants of tt and K through third order in this 
expansion compared to the model exact results obtained 
earlier via direct analytic continuation of the DSE so- 
lutions for propagators. Agreement to three significant 
figures is obtained at second order for tt while there is ev- 
idently still an error of about 2% in fx after third order 
for K. 

For the p, we show in Fig. || the BSE eigenvalue A(P 2 ) 
for the first three orders of the Taylor expansion. We 
also display the model-exact behavior obtained by direct 
integration throughout the relevant complex domain. At 
the physical point A = 1, a converged value of the mass is 
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evident in the figure at second order in the Taylor expan- 
sion method. However, the figure also suggests that at 
larger timelike P 2 a converged behavior is not yet clearly 
evident at third order. 
The results in Table 



III 



for p, K* and <j) illustrate this 
in more detail. At second order the masses and f p are 
within 1% of the model exact results while fx* and are 
within 5% and 3% respectively. Since K* is not a charge 
conjugation eigenstate, the solution of the BSE is more 
difficult due the the lack of definite parity in the variable 
q ■ P. Furthermore, one expects larger errors in the de- 
cay constants than in the masses since the former require 
normalization of the BS amplitudes via Eq. (P) and the 
derivatives therein add to the numerical difficulty. The 
third order results for the masses improve slightly for p 
and K* but m</, deteriorates slightly to acquire a 2% er- 
ror; the decay constants remain within 5% of the exact 
values. The results do not indicate that the Taylor ex- 
pansion has yet converged and this correlates with the 
behavior evident in Fig. 0. Table III also makes clear 



that a truncation of the eight Dirac covariants to just 
the dominant one (7^) will lead to a persistent error of 
at least 15-20%. In the case of the K* this tends to raise 
the mass to the level where we are unable to determine 
a converged result. 




-0.5 



0.5 



P (GeV ) 



FIG. 3: The BSE eigenvalue in the p channel calculated at 
the indicated orders in the Taylor expansion compared to the 
model exact result ||. The physical point is at A = 1. 

The evident difficulty in use of the Taylor expansion 
approach as the mass of the state increases is most 
likely related to the occurrence of complex conjugate 
singularities of the employed quark propagator ampli- 
tudes a a (q±). In rainbow approximation, which we em- 
ploy, the occurrence of complex conjugate logarithmic 
branch points with timelike values of is well-known 



for the fcrmion propagator in both QED [g4| 25 and 
QCD [^6[ |?J. For the present model, we find numeri- 
cally that the non-analytic behavior nearest to the origin 
is at q 2 ± = -0.207 ± zO. 331 GeV 2 for the u/d quark, and 
at q\ = -0.376 ± i0.602 GeV 2 for the s quark. These 
points are outside the integration domains needed for 



TABLE III: The vector meson masses and electroweak decay 
constants (in GeV) . Results are shown through several orders 
in the Taylor expansion treatment of q\ in quark propagator 
amplitudes. The effect of truncation of the 8 covariants to just 
the canonical (Dirac matrix) covariant is indicated. Compar- 
ison is made with results from use of quark amplitudes along 
5f(5±) only, and with the model-exact results. 





P 




K* 




4> 






m p 


f P 




fK* 




U 


1 st order Taylor 














All 8 amplitudes 


.669 


.222 


.848 


.276 


1.009 


.299 


Vi only 


.811 


.255 






1.202 


.331 


2 nd order Taylor 














All 8 amplitudes 


.732 


.209 


.949 


.253 


1.064 


.266 


Vi only 


.856 


.214 






1.235 


.271 


3 rd order Taylor 














All 8 amplitudes 


.735 


.196 


.934 


.253 


1.054 


.253 


Vi only 


.844 


.187 






1.182 


.213 


Real axis only 


.711 


.217 


.872 


.241 


1.066 


.288 


Model exact [9 


.742 


.207 


.936 


.241 


1.072 


.259 


Experiment 


.770 


.216 


.892 


.225 


1.020 


.237 



the model-exact p and <j) states respectively. However 
the critical bound state masses for which these singu- 
lar points of the propagators would just begin to enter 
the integration domain are 1.09 GeV for a uu state, and 
1.47 GeV for a ss state. A Taylor expansion should di- 
verge at these points; it is probably the precursor of this 
that is beginning to show for the p in Fig. |3|. 



C. Strong Decays 

Since the p and <\> appear as resonance poles in the 
timelike behavior of the charge form factors of the ir and 
K mesons, the pole residues involve the coupling con- 
stants g p7 nv and g^KK- If the charge form factors, or 
more generally the jtttt and jKK vertex functions, are 
formulated in impulse approximation, then the values of 
<7 p7r7r and g^KK extracted from the pole residues will be in 
impulse approximation. The quality of the latter values 
will depend directly upon the quality of the strengths 
of the 77T7T and "fKK vertex functions in the relevant 
timelike region. It is known that the impulse approxima- 
tion for electromagnetic coupling to mesons conserves the 
electromagnetic current (and provides the correct charge) 
independent of model parameters as long as the meson 
BS amplitudes and photon-quark vertex are in ladder 
approximation and the quark propagators are in rain- 
bow approximation fi"o|| . We assume that the dynamics 
which produces the correct strength at the photon point 
will also produce a good quality strength at the mass 
shell of the vector mesons and hence render the impulse 
approximation good for g p7rK and g^KK- The extension 
of this argument to cover the decay K* + — > K + ir° is af- 
forded by the study |ll[] of the semileptonic K(13) decay 
K + — > lvi"K° where K* + appears as a pole in the W + 
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TABLE IV: Vector decay coupling constants calculated with 
the imaginary parts of amplitudes treated through order two 
in the Taylor expansion. The dependence on the number of 
invariant amplitudes employed is indicated. The K* decay 
process shown in brackets is related by isospin symmetry to 
the former process by a factor 1 / \j2. The notation (v,p) indi- 
cates the number of invariant amplitudes used for the vector 
meson and for the pseudoscalar mesons respectively. 





(v,p)=(8,4) 


(v,p)=(5,4) 


(v,p)=(l,l) 


Expt 


Qp — >7T7V 


5.14 


4.93 


8.8 


6.02 


g^KK 


4.25 


4.06 


6.91 


4.64 




4.81 


4.56 




4.60 


(fe* + ^iC+7T°) 


(3.40) 


(3.22) 




(3.26) 



vertex. 

Here the coupling constants for the decays of p, K* 
and 4> to a pair of pseudoscalar mesons are calculated in 
impulse approximation according to Eqs. (0-^). The con- 
straints on the two external Euclidean momenta P and 
Q needed to satisfy the three mass-shell conditions entail 
an analytic continuation to a complex value for at least 
one component. For example, in the p — > Ttir case, with 
the momenta defined in Eq. (jl]), we have Q = (im p ,0), 
P 2 = m 2 /A~ m\ and P ■ Q = 0. A by-product of this 
is that two of the quark propagators and two of the me- 
son BS amplitudes are needed in domains of the complex 
q 2 -plane. In all cases we employ second order Taylor 
expansions about the closest point on the positive real 
axis. We also test the dependence upon the number of 
Dirac covariants retained for the vector and pseudoscalar 
BS amplitudes. With a maximum of eight for the vector 
and four each for the pseudoscalar BS amplitudes, there 
are a maximum 128 distinct quark loop integrals to be 
performed. 

The results are summarized in Table |y| in comparison 
with experimental values obtained from decay widths. 
With just the dominant covariant employed for each me- 
son (75 f or p seudoscalars, 7^ for vectors, i.e, (v,p)= (1,1) 
in Table pv| ) the coupling constants that can be obtained 
are about 50% larger than experiment. It is known from 
work on the vector meson BSE that a very efficient trun- 
cation for soft physics can be obtained with a particular 
set of 5 covariants ||. With use of this truncation to- 
gether with all 4 covariants for each pseudoscalar, the 
results in Table IV indicate that the dominant physics 

is 



has been captured; the K* decay is within 1%, the 
within 12% and the p decay is 18% less than experiment. 
The addition of the remaining three transverse vector co- 
variants leads to modest improvement giving a deviation 
from experiment of between 5% and 10% with the error 
being larger if the vector meson is lighter. 

Since the width of the p is almost 20% of its mass 
while the widths of the 4> and K* are significantly less 
important, we expect the ladder approximation for the 
BSE kernel (which omits the strong channels 7T7t, KK 
and Kir respectively) to be less accurate for the p than 
for the 4> and K*. Accordingly we speculate that this is 



TABLE V: The present results for the coupling constants 
g v ^p P (calculated with all covariants) compared to values ex- 
tracted from a pole fit to timelike electroweak form factors, 
and also to experiment. 



Qv^pp 


this work 


pole fit p8|, p9| 


Expt 


Qp — -'7T7V 


5.14 


5.2 


6.02 


94>—*KK 


4.25 


4.3 


4.64 




4.81 


4.1 


4.60 



largely the reason why the result for g p -^ V7r in Table [TV 
deviates from experiment twice as much (15%) as do the 
other decay constants. 



The results shown in Table IV compare well with 
coupling constants extracted from the behavior of elec- 
troweak form factors calculated, independently and in the 
same model, at timelike momentum near the vector me- 
son poles |28|, ^9] . The general behavior of a pseudoscalar 
meson charge form factor Fp(Q 2 ) in this domain is p8| 



F P (Q 2 ) 



9v- 



9v {Q 2 + ml - im v f„) 



(20) 



Here m 2 /g v is the u-7 coupling strength calculated from 
the D->e + e~ decay, and T v is the vector meson width 
associated with the mass-pole. In the present work the 
ladder approximation to the BSE produces real mass- 
poles. The vector decay coupling constants obtained in 
this way from a fit to the calculated F 7T (Q 2 ) and Fk(Q 2 ) 
form factors and from the form factor for the semilet 
tonic K(13) decay K + — > Iviir are shown in Table 
These results from a pole fit employ all relevant covari- 
ants for BS calculations and employ complex plane con- 
tinuations where necessary. From this perspective the 
comparison with the left-most column of Table IV pro- 



vides a rough measure of the effectiveness of the Taylor 
expansion method for BSE solutions. 



V. ESTIMATE OF HEAVIER STATES 

In the present formulation, the nearby singularities dis- 
cussed in Sec. IV B hinder a direct application of our 



rainbow-ladder model for masses above about 1.2 GeV. 
In particular, it is of interest to go beyond the present 
S-wave states to consider orbital excitations such as the 
axial vectors a\ and b\ . One can however provide an indi- 
rect estimate as follows. Although the quark propagator 
amplitudes <r a (q±) are complex- valued in the complex q\ 
plane, the ladder BSE eigenvalue A(P 2 ) is real and so is 
the mass. The evident cancellations may allow the real 
axis behavior of the quark amplitudes to dominate the 
outcome. As a test, we introduce a real-axis approxi- 
mation by setting the imaginary part of the momentum 
argument of the quark propagator amplitudes a^q 2 .) to 
zero in the BSE kernel. Then, instead of Eq. dl9|), one 
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has 



M±) « a a {q 2 -m 2 /A) , 



(21) 



where a = v,s and q 2 > 0. This is not quite identical to 
zeroth order in the Taylor expansion of Eq. ([l9]) because 
the latter expands about the nearest point on the posi- 
tive real axis, while in Eq. ( pl| ) the interval [— m 2 /4, oo] 
is used. Thus the required domain of the timelike real 
axis is treated exactly and this can be important for the 
heavier mesons. 

To implement this approximation, we use the quark 
DSE to continue the propagator solutions to the timclike 
(negative real) p 2 axis. In this respect, it is advantageous 
to enhance numerical stability of the present model by 
modifying the infrared behavior of the factor T(k 2 ) in 
the small second term of the effective coupling, given in 
Eq. (|l8|), so that, like the first term, it is identically zero 
at k 2 — 0. Its UV limit and range should remain much 
the same. To achieve this we use 



1 _ e -(«/4l 



(22) 



With the real axis approximation of Eq. the result- 
ing modified model yields mass estimates that tend to be- 
come more accurate as the mass increases. From Table 
m-x/K are underestimated by 10% and fn/K are under- 
estimated by 20% compared to the model-exact values. 
However, as shown in Table HI, the produced ground 



state vector meson masses m p , ttik* and are within 
4%, 7% and 1% of the model-exact values respectively; 
the decay constants deviate by 5%, 0% and 11% respec- 
tively. We therefore expect this real-axis approximation 
to provide mass estimates in the 1-2 GeV range with an 
error of about 5-10%. 

For a u/d quark axial vector state, the required general 
form of BS amplitude is simply 75 times Eq. (|l^). For 
the ai (1 ++ ), amplitudes V5, V7 and Vg will be odd in 
q ■ P and the rest will be even; for the opposite C-parity 
state b\ (l -1 ) the opposite is true. The real axis approx- 
imation gives the following axial vector estimates: m ai ~ 
0.891 GeV and m bl sa 0.775 GeV. Hence the orbital ex- 
citation energy m ai — m p is 0.150 GeV and thus a factor 
of 3 too small in this model. Indications from estimates 
attempted using the complex plane information are sim- 
ilar, as are recent results from a related study [Q. Mod- 
els of the present rainbow-ladder type are significantly 
too attractive for these orbital excitations. Other studies 
of a\ and b\ based on the DSEs have used a separable 
approximation where the quark propagators are the phe- 
nomenological instruments |32|, |3^| and these studies 
find more acceptable masses in the vicinity of 1.3 GeV. 
The manner in which the effective interaction is modeled 
is clearly important. Preliminary results from an exten- 
sion of the present rainbow-ladder level through a 1-loop 
dressing of the quark-gluon vertex, while preserving the 
vector and axial vector WTIs, indicate that the sepa- 
ration of the axial vector states from the vector states 
increases significantly and becomes quite acceptable [p4j . 



Recently there has been interest in meson states with 
quantum numbers that are called exotic in the sense 
that they cannot be produced as qq bound states within 
static quantum mechanics where C-parity is given by 
C = {-l) L+s . The tti(1400) state, formerly called the 
p(1405), is a resonance with J c = 1 h that is sup- 
ported by evidence from nd scattering and pd annihila- 
tion ]lq] . In a covariant field theory treatment of bound 
states via the BSE, the extra degree of freedom repre- 
sented by relative time, or relative energy, allows states 
to have an additional classification in terms of an asso- 
ciated time-parity quantum number k and C-parity is 
given by the more flexible form C = ( — 1) L+S+K . The 
states with k = odd have no static quantum mechanical 
qq analog and in this sense the exotic J PC obtainable 
from the ladder BSE are evidently a realization of the 
multi-particle features of covariance. Although a full un- 
derstanding within field theory of the norm of these odd 
solutions of the BSE remains to be achieved J3^|, ex- 
amples from otherwise realistic models are of interest in 
connection within on-going experimental searches. 

A covariant separable BSE model, closely connected 
with the present model, has recently produced |3j| 
m 7ri = 1.439 GeV, i.e., about 100 MeV above the 01, 
which in that model is m ai = 1.337 GeV. To estimate 
whether this is consistent with the present study, we use 
the real axis approximation to search for a J PC = 1 
solution. The general form for such a BS amplitude is 
the same as Eq. ( |l5|) for the 1" amplitude except that 
one reverses the odd-even property of the invariant am- 
plitudes in q ■ P. That is, V3 and Vq are to be even, the 
rest are to be odd. This produces a tti solution 112 MeV 
above the a\ within the same approach. In this respect 
we agree with Ref. Q. We speculate that if the prob- 
lem of the ground state orbital excitations being about 
300 MeV too low is remedied, then the tti state in this 
model should rise by a similar amount to be around 1.3- 
1.4 GeV. 



VI. 



DISCUSSION 



We have studied vector meson strong decays within 
a Euclidean space model of QCD based on the Dyson- 
Schwinger equations truncated to ladder-rainbow level. 
The infrared structure of the ladder-rainbow kernel is 
described by two parameters; the ultraviolet behavior 
is fixed by the one-loop renormalization group behav- 
ior of QCD. Within the u/d and s quark sector we have 
obtained the coupling constants for the strong decays: 
p — > 7T7T, K* — > Kir, and (j> — > KK. The deviation from 
experiment is between 5% and 10% with the error being 
larger if the vector meson is lighter. 

Our method of calculation employed a Taylor expan- 
sion to continue the dressed quark propagators from the 
real spacelike p 2 axis into the domain of the complex 
plane required by the mass-shell condition of the meson 
bound states. We used this both for the solution of the 
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BSE and for calculation of decays in the impulse approx- 
imation. The benefit is that knowledge of the dressed 
quark propagators is needed only along the real space- 
like p 2 axis and this simplifies the numerical requirements 
considerably. The vector and pseudoscalar meson masses 
compared well with the model-exact results obtained pre- 
viously by direct continuation of the dressed quark prop- 
agators into the complex plane via the DSE. In particu- 
lar, the masses and electroweak decay constants at sec- 
ond order are all within 5% of the model-exact results. 
The increased deviation observed at third order and for 
heavier masses is attributed to nearby complex conjugate 
singularities in the propagators that are known to occur 
in rainbow solutions of the DSE. 

For mesons heavier than the <fi, these singularities hin- 
der the applicability of Euclidean space implementations 
of the DSEs as presently formulated. The Taylor ex- 
pansion method cannot lessen this difficulty. We there- 
fore explored a real-axis approximation that exploited the 
cancellations of imaginary parts evident in the way the 
complex- valued quark amplitudes off the real axis enter 
the BSE kernel to produce real eigenvalues and masses. 
The real-axis approximation avoids complex plane sin- 
gularities and, based on its performance for the pseu- 
doscalar and vector mesons, we conclude it can be used 
for mass estimates in the 1-2 GeV region with an ac- 
curacy of about 5-10%. We then produced estimates of 
m ai and for orbital excitations as well as m 7ri for 
the proposed exotic 1 h vector state. We mentioned ev- 



idence that an extension beyond ladder-rainbow level to 
include dressing of the quark-gluon vertex could address 
our finding that the present model is too attractive for 
states just above the ground state vectors. With this, 
our estimate for m 7Tl is 1.3-1.4 GeV. 

Work in progress |35) on the case of propagators hav- 
ing explicit complex conjugate poles, indicates that a re- 
formulation produces well-defined integrals and integral 
equations above the critical masses. Beyond rainbow- 
ladder approximation, very little is known about the sin- 
gularity structure. In one simplified model that includes 
a fully dressed quark-gluon vertex in a way that has al- 
lowed exploration of analytic properties, the non-analytic 
behavior is relegated to essential singularities at infin- 

ity §• 
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